Draft version March 8, 2013 

Preprint typeset using L£-T^X style cmulateapj v. 5/2/11 



(N 

o 

(N 
C 

(N 



< 
O 

Oh: 
6 

CO 

o3 



> 

O 

(N 



13 



AN UPPER LIMIT TO THE VELOCITY DISPERSION OF RELAXED STELLAR SYSTEMS WITHOUT 

MASSIVE BLACK HOLES 

M. Coleman Miller 1 and Melvyn B. Davies 2 
Draft version March 8, 2013 

ABSTRACT 

Massive black holes have been discovered in all closely examined galaxies with high velocity dis- 
persion. The case is not as clear for lower-dispersion systems such as low-mass galaxies and globular 
clusters. Here we suggest that above a critical velocity dispersion ~ 40 km s^ 1 , massive central black 
holes will form in relaxed stellar systems at any cosmic epoch. This is because above this dispersion 
primordial binaries cannot support the system against deep core collapse. If, as previous simulations 
show, the black holes formed in the cluster settle to produce a dense subcluster, then given the ex- 
tremely high densities reached during core collapse the holes will merge with each other. For low 
velocity dispersions and hence low cluster escape speeds, mergers will typically kick out all or all 
but one of the holes due to three-body kicks or the asymmetric emission of gravitational radiation. 
If one hole remains, it will tidally disrupt stars at a high rate. If none remain, one is formed after 
runaway collisions between stars, then it tidally disrupts stars at a high rate. The accretion rate after 
disruption is many orders of magnitude above Eddington. If, as several studies suggest, the hole can 
accept matter at that rate because the generated radiation is trapped and advected, then it will grow 
quickly and form a massive central black hole. 

Subject headings: Accretion, accretion disks — Black hole physics — Galaxies: clusters: general — 
Galaxies: bulges — Gravitation — (Stars:) binaries: general 



1. INTRODUCTION 

Observations over the last two decades have revealed 
central massive black holes in all sufficiently wel l- 
observed massive galaxies (e.g., iGultekin et al.l 1201 ID . 
However, the case is not as clear for lower-mass galax- 
ies or globular clusters, and indeed although there 
is evidence for black holes in some low-mass galaxies 
(|Greene et al.l 120101 : IKuo et ail 120 111) there are exam- 
ples of galaxies that clearly do not have black holes 
that follow the st andard mass - ve l ocity dispersion 
(M - a) relation (|Merritt et al.l 120011 : iGebhardt et al.1 
2001) and the case for globular clusters is far from 
clear (e.g., iGerssen et all l2002t iMcNamara etafl 120031: 
iBaumgardt et al.ll2003l: iStrader et al.ll2012fl . 

Here we approach this question by focusing on the ve- 
locity dispersion rather than the mass of a stellar sys- 
tem. In Section 2 we show that above a critical velocity 
dispersion <7 cr ;t ~ 40 km s _1 , the total binding energy 
in primordial binaries that can be tapped in three- and 
four-body interactions is significantly less than the total 
binding energy of the system as a whole, and hence if such 
systems are dynamically relaxed they will undergo deep 
core collapse essentially unhindered by dynamical heat- 
ing from b inaries (thus leading to o ne of the scenarios dis- 
cussed bv lBegelman fc R ces 1978jin the context of more 
massive clusters). We note that the galaxies seen thus 
far without massive black holes have velocity dispersions 
below this limit (e.g., NGC 205 has a = 39 km s _1 and 
M33 has a = 24 km s _1 ; see lGultekin et al. 1120091 and ref- 
erences therein). In Section 3 we discuss the evolution of 
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binary-free systems. Previous studies have demonstrated 
that the black holes in such systems sink rapidly to the 
center and interact mostly with each other in a dense sub- 
cluster. This leads to three paths, all of which culminate 
in the formation of a massive black hole: (1) For suffi- 
ciently high escape speed systems dynamical interactions 
result in runaway merging of the black holes into a mas- 
sive hole. For lower escape speed systems either one or 
zero black holes remain after ejection of merged pairs due 
to asymmetric emission of gravitational radiation dur- 
ing coalescence or Newtonian recoil from interactions of 
black holes with dynamically formed binaries. (2) If one 
black hole remains then it tidally disrupts ordinary stars 
and consumes the remnant disks quickly, hence grows 
rapidly into a massive black hole; other growth mecha- 
nisms, such as the accretion of nascent gas or winds, are 
insignificant. (3) If no black holes remain then runaway 
collisions form a massive star that evolves into a black 
hole, and this first black hole grows via accumulation of 
tidally disrupted stars. Thus once binary support is re- 
moved, massive black hole formation is assured as long 
as holes consume tidal remnants quickly. In Section 4 we 
determine the minimum mass of a black hole formed via 
these paths and discuss the implications of this scenario. 



2. VELOCITY DISPERSION THRESHOLD FOR DEEP 
CORE COLLAPSE 

Stellar systems that are in virial equilibrium evolve via 
two-body interactions over their relaxation time, which 
for a star of mass m in a system of velocity dispersion a 
at a location with an average stellar mass density p is 
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(| Spitzerl 1 1 987h . where In A - 5 - 10 is the Coulomb log- 
arithm. The evolution of an isolated stellar system is 
towards a greater concentration of stars in the center 
balanced by a greater expansion of the cluster on the 
outskirts; there is a productive analogy with thermody- 
namics, in which this behavior can be seen as the grad- 
ual increase of cluster entropy (the greater phase space 
accessed by the outer stars more than makes up for the 
diminished phase space accessed by the stars in the core) . 
It was demonstrated several decades ago that if all the 
stars are single (as opposed to being in binary or multi- 
ple systems), then over a timescale that scales with the 
relaxation time at the half-mass radius for a typical star 
(where the multiple is ~ 15 for an initially Plummer 
sphere of equal-mass stars but is ~ 0.2 if there is a broad 
initia l mass function; see IPortegies Zwart fc McMillan 
2002), the core becomes so dense that it loses thermal 
contact with the rest of the cluster and the core under- 
goes a collapse such that the nu mber density in the inner 
portions scales a s n ~ r~ 22 (|Lvnden-Bell fc Eggletonl 
119801: iCohnl 11980ft . If we take present-day nucl ear star 
cluste rs as an example, then from Figure 1 of iMerrittl 
(2009) we find that most have half-mass relaxation times 
less than few x 10 10 yr and thus are candidates to col- 
lapse within a Hubble time if they had broad initial mass 
functions and no central massive object to supply energy. 

Binaries are the key to sustaining a cluster against this 
collapse. When number densities become high enough 
that binary-single interactions are common, such interac- 
tions can harden the binary and hence inject energy into 
the cluste r that decreases its densit y. Many calculations 
(see, e.g.. lHenonlll961l : lHeggidll975l for pioneering work) 
have shown that binaries that are initially hard (meaning 
that their binding energy exceeds the kinetic energy of a 
typical single star) tend to harden via binary-single inter- 
actions, whereas initially soft binaries tend to soften and 
eventually break up. Consistent with this expectation, 
globular clusters have a significant ly smaller binary frac- 
tion t han is seen in the field (e.g.. (Rubenstein fc Bailynl 
ll997tlMilone et al.ll2012h . 

In principle, even a very small number of binaries could 
have enough binding energy to hold off the collapse of a 
cluster. Consider for example a reasonably rich globu- 
lar cluster with a velocity dispersion of 10 km s _1 . A 
binary of two solar-mass stars near contact, with an or- 
bital radius of 0.01 AU, has ~ 10 3 times the binding 
energy per mass that a single cluster star has in kinetic 
energy, so if 0.1% of stars are in such binaries the energy 
to hold off cluster collapse appears to be present. White 
dwarfs are 100 times smaller yet, so it might seem that 
if there is one near-contact white dwarf binary in a clus- 
ter of 10 5 stars its binary interactions could successfully 
oppose core collapse. 

This is of course not true, for two reasons. First, as the 
scmimajor axis of a binary shrinks, its close interactions 
with single stars have a greater and greater chance of 
destroying the single star or one of the binary stars, hence 
the kinetic energy o f recoil is not shared with the cluster 
(|Davies et al.lll~9 94). As an example, a tight white dwarf 
binary cannot eject a main sequence star in this way. 
Second, even if a three-body interaction is clean, a star 
that is thrown completely from the cluster cannot share 
its kinetic energy with the cluster and the only expansion 



of the core comes from the comparatively minor effect 
that the core now has lost one star's worth of mass. 

The available binding energy from binaries is thus 
limited; clusters having higher velocity dispersions hav- 
ing a more limited available binding energy. As we 
now argue, this means that above a velocity dispersion 
o"crit ~ 40 km s , the binaries cannot hold off core col- 
lapse. It should be noted that the velocity dispersion of 
a cluster will evolve as a function of time, with velocity 
dispersions being somewhat larg er in the past when the 
cluster was more massive (e.g., iGiersz fc Heggiel [200l 
iKiipper et al.l 120101 ) . The effect could be particularly 
enhanced for clusters containing multiple stellar popu- 
lations where a large fraction of t he first generation of 
stars are lost (|D'Ercole et al.ll2008l ). However the veloc- 
ity dispersion at later times will be more relevant to the 
discussion in this paper, as this is when core collapse may 
typically be possible (i.e., on timescales longer than the 
half-mass relaxation time). 

As a first estimate of the available binding energy for a 
binary with initial semimajor axis ao, we assume that the 
eccentricity distribution of binaries with a given semi- 
major axis is a thermal distribution P{e < eo) = e 2 , 
truncated at the maximum eccentricity e max allowed for 
pericenter distances greater than some minimum r p . m ; n 
(this could be the pericenter distance at which stars col- 
lide), a(l — e max ) = r p m i n for a semimajor axis a. Thus 
a fraction of orbits are allowed, hence the binding 
energy that can be released from semimajor axis a + da 
to a is weighted by e^Ja) = 1 - 2r p , min /a+ (r Pimin /a) 2 . 
Thus the total available binding energy from an initial 
semimajor axis ao with stars of mass m is 
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so for simplicity we will approximate i?bind,tot as zero be- 
low 10r Pjm i n and Gm 2 /6r Piin i n above it. 

Our next step is to note that for stars formed in a 
low-density environment, there is roughly one binary 
per single star, and the binary semimajor axes are ap- 
proximatel y equally distribute d in In a from 0.01 AU to 
~ 10 4 AU (Popova et al. 1983). In an environment where 
binaries beyond a certain semimajor axis are ionized by 
binary-single encounters, the fraction of binaries will be 
decreased. For example, if we begin with six single stars 
and six binaries and ionize the ones larger than 1 AU, 
we now have fourteen single stars and two binaries. If 
as above we now only concentrate on the binaries larger 



than 0.1 AU = 10r D 



this represents / ~ 7% of the 



stars in the system. Thus the binary binding energy per 
all stars in the system is 

e bin /star = fGm 2 /(6r Ptluin ) . (4) 

This is to be compared with the binding energy per star 
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in the cluster, which by the virial theorem equals the 
kinetic energy per star in the cluster, or 

edustcr/star = \^na 2 (5) 

for a velocity dispersion a. The point at which 
ebin/star < e c i ustcr /star is the point at which core col- 
lapse is theoretically possible. From the numbers above, 
if the single stars and the binary components both have 
masses w lAf Q this happens when a ~ 40 km s~ x , 
meaning that interactions with binaries of semimajor 
axis <; 0.5 AU have positive total energy and are thus 
soft, core collapse can proceed. If the initial distribution 
of binary binding energies is extremely unusual, e.g., if 
most stars are formed in binaries with semimajor axes 
less than 0.5 AU, then the supply of binary energy would 
be greater and the threshold velocity dispersion could in 
principle be raised. Barring such an unexpected distri- 
bution, however, th e threshold should b e robu st. 

Indeed, work by iChernoff fc Huang] (JT996) suggests 
that there may be less binary energy available than we 
derive above. They take into account that, rather than 
simply resetting the eccentricity of a binary, a binary- 
single encounter can be resonant and hence for a given 
interaction there is a greater chance to get to a very 
small separation. From their Figure 4 we infer that for 
solar- type stars and a = 40 km s _1 a typical energy 
AE « 6 x 10 46 erg can be extracted from an initially 
hard binary, whereas equation ([3]) gives roughly an order 
of magnitude larger energy. Thus at a — 40 km s _1 , 
and perhaps at a slightly lower velocity dispersion, the 
energy that can be extracted from primordial binaries is 
significantly less than the binding energy of the cluster, 
hence such clusters can undergo core collapse without 
being impeded significantly (three-body binary forma- 
tio n and two-body tida l ca pture are also insig nificant; 
see iHut fc Bahcalllll983l and lFabian et alJll975L respec- 
tively). 

3. PATHS TOWARDS MASSIVE BLACK HOLE FORMATION 

We now evaluate the paths towards massive black hole 
formation that we mentioned in the introduction: run- 
away merging of black holes, tidal disruption of stars 
by a single remaining black hole, and formation of a 
new black hole from runaway collisions of stars, fol- 
lowed by tidal disruption of stars by that black hole. 
The hrst important question is one of time scales. In 
all three paths, the overwhelmingly longest phase is the 
initial progression to core collapse. To see this, note 
that the time to core collapse is a multiple less than 
unity (~ 0.2 for systems with a broad ma ss distribu- 
tion; see IPortegies Zwart fc McMillan! 120021 ) of the re- 
laxation ti me of the nuclea r star cluster, which from Fig- 
ure 1(a) of lMerrittl ([27)09h is i rbt ~ 10 9 yr(Af/10 6 M ) 
with a spread of a factor ~ 10. There is some evidence 
that nuclear star clusters obey a similar M — a rela- 
tion to that seen for higher-mass black holes. There is 
some observational e vidence for this; e.g., Figure 2 from 
iFerrarese et al.ll2006( ) indicates that nuclear star clusters 
might have the same M — a slope as has been found for 
black holes (see IGultekin et ahl 120091 for a recent discus- 
sion of this relation) but offset so that the mass is a factor 
of ~ 10 higher than the black hole mass would be. If we 
loosely equate the velocity dispersion of the cluster with 



that of the surrounding bulge, this gives a cluster mass 
of 

M d » 10 6 M Q (er/40 km s^ 1 ) 4 (6) 

based on the scalings of IGultekin et al.l (|2009l ) . Thus 
clusters with a Sa 100 km s~ have a chance to undergo 
core collapse within a Hubble time. Figure Q] illustrates 
the paths we consider. 
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Fig. 1. — Paths towards massive black hole formation in a stellar 
cluster. At velocity dispersions less than ~ 40 km s — 1 , heating 
from binaries prevents full core collapse. At velocity dispersions 
greater than ~ 100 km s^ 1 , a typical nuclear star cluster will have 
too long a relaxation time for its core to collapse in a Hubble 
time, although a massive black hole can form in other ways. At 
intermediate velocity dispersions, full core collapse will occur and 
will likely result in either zero or one remaining stellar-mass black 
hole. In the latter case, the hole will grow via tidal disruption of 
stars; in the former, the stars will undergo runaway collisions that 
produce a black hole, which will then grow via tidal disruption of 
stars. 

This step is necessary for all three paths we discuss. 
For a collapsed core, t h e self -similarity arguments of 
Lynden-Bell & Eg gletonl (|1980| ). and the classic simula- 
tions of lCohrd (|1980E show that if all the stars are treated 
as point masses and no three-body binary formation is 
allowed, then the density of a single-mass system evolves 
towards a n oc r~ 2 2 configuration. This is quite close to 
a singular isothermal sphere n cx r~ 2 , hence we will as- 
sume that the velocity dispersion is nearly constant in the 
collapsed region. As a result, the relaxation time scales 
roughly as so the evolution timescale is shorter by 
orders of magnitude in the core of the cluster than it is 
in the cluster as a whole. 

As a result, once the core collapses all three paths 
are traveled in a time much shorter than the time to 
collapse. For example, runaway mergers between black 
holes (or, in the third path, runaway collisions between 
stars) occur roughly on the core relaxation timescale, be- 
cause when the number density is not yet sufficient for 
frequent mergers or collisions, further relaxation will in- 
crease the density on the relaxation time until interac- 
tions are frequent. Thus the only limiting factor is the 
initial collapse time. We also note that unlike in the sce- 
nario of r unaway collapse of young massive cl usters pro- 
posed by IPortegies Zwart fc McMillan! ([2002) , the time 
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window for runaway collisions of stars to form a single 
black hole (in the third path) is not millions of years, 
but billions of years. The reason is that when the cluster 
is young enough that initially all stars are on the main 
sequence, supernovae from the most massive stars begin 
at ^2.5 Myr and proceed for many stars, causing the 
core to lose a large amount of mass to the ejecta and 
therefore expand and lower the number density. In con- 
trast, in our picture the evolution to core collapse is much 
later, perhaps billions of years, hence the remaining stars 
are low-mass and thus only the collision product will be 
massive enough to explode; very little mass is lost, so the 
density remains high. 

In addition to the general core collapse, in a multi- 
mass system there is considerable mass segregation. This 
means that the stars in the core will tend to be to- 
wards the massive end, perhaps ~ 1 M after billions 
of years. In addition, of the objects likely to be present 
after a long time, stellar-mass black holes will be by a 
factor of a few to several the most massive. Many stud- 
ies (e.g., iMackev et al.l I2008T) have concluded that the 
black holes then form a dense subcluster in which the 
holes interact mainly with themselves. If, as in our sce- 
nario, there are no binaries, then the holes can reach ex- 
tremely high density in the center of the subcluster and 
capture each other via emission of gravitational radia- 
tion in initially hyperboli c two-body encounters. From 
IQuinlan fc Shapfra ([1989D . the critical pericenter for a 
two-body gravitational wave capture between two black 
holes with a total mass M = m\ + vni and a reduced 
mass \i = mirriz/M is 

r p . GW = 8.5 x 10 8 cm (M/20 M ) 5 / 7 ( 
x(/V5 M ) 2 / 7 (a/4O km s^ 1 )- 4 / 7 , [ ' 

and hence their gravitationally focused cross section is 

E bh = 2nr p GM/a 2 « 9 x 10 23 cm 2 (M/20 Mq) 12 ' 7 
x(ji/5 M Q ) 2 / 7 (cr/40 km s" 1 )- 18 / 7 . 

(8) 

When two black holes capture each other in this way, 
their inspiral is extremely rapid: from [Peters (1964]), the 
inspiral time is 

5 4/t 2\7/2 

if j I j.\ 5 c a (1 — e ) ' 

a/yaa/ai) — 64 G3flM 2 (i + 73 e 2/ 2 4+37e 4 /96) (9) 
« 10 5 yr(a/i? ) 4 (l - e 2 ) 7 / 2 

where the approximation is for e ~ 1 and our given num- 
ber assumes mi = m 2 = 10 M©. For a « li? 
7 x 10 10 cm and e > 0.99 (so that r p < r PiG w), the 
inspiral time is therefore less than 0.1 years, and for a 
fixed r p the inspiral time scales as a 1 / 2 so that even for 
an initial a = 100 AU the inspiral time is just a few years. 

When the holes do merge they emit gravitational ra- 
diation that is in general asymmetric, meaning that 
the remnant single black hole will recoil relative to its 
original center of mass. Studies of black hole recoil 
(iBaker et al.l2007l[2 008; Lo usto fc Zlochowerl2008ll200l 
Lousto et al.l I2010L Ivan Meter et al.l 120101 : lLousto et al.l 
20121 ) show that although kicks from the coalescence of 
nonspinning black holes are limited to < 200 km s , 
rapidly spinning black holes can produce remnants that 
travel at thousands of kilometers per second relative to 
their original center of mass. Thus in this environment, 



unlike in the conditions th at may exist in the z > 10 
universe (|Davies et al.ll20rl . mergers that are restricted 
to comparable-mass black holes are most likely to lead 
to an ejection of the remnant. 

As pointed out to us by S. Sigurdsson (private com- 
munication), for low velocity dispersions the ejection 
of black holes is likely to be dominated by encounters 
with hard binaries formed by the interact ion of three ini- 
tially hyperbolic black holes. iHutl (|1985l ) finds that the 
rate of formation of "immortal" binaries by this process 
(i.e., binaries that are not later softened and ionized) is 
TI3B = 126G 5 m 5 n 3 /<7 9 , where he assumes objects of iden- 
tical mass m and a is the three-dimensional velocity dis- 
persion. Thus the ratio of the formation rate per volume 
of these binaries to the rate of gravitational wave capture 
of black holes by each other, assuming equal masses, is 

-7———— — « 200(n/10 10 pc- 3 )(m/10 M ) 3 (cr/4O km s" 1 )" 52 / 9 . 

"■BHBH 

(10) 

Thus for low to moderate velocity dispersions, and high 
number densities, binary-single ejections are likely to 
dominate. The result will be similar to the case in which 
only double black hole mergers occur: there will either 
be zero or one hole left. For the rest of this section we 
concentrate on ejections by mergers. 

We set up a simple simulation of the evolution of a 
black hole subcluster with no binaries. We assume that 
initially there are either 100 or 101 black holes; note 
that even if all mergers eject the remnant, having an 
odd number initially guarantees that one will survive be- 
cause with no binaries the interactions are pairwise. The 
distribution of black hole masses is not well-established, 
and the distribution of their spins is even less so, but 
as an illustrative example we draw the initial masses of 
the black holes from the range [5, 30) M , with a distri- 
bution dN/dM oc M -2 , and the initial spins are drawn 
uniformly from the range cJ/(GM 2 ) = [0, 1). 

We simulate the evolution of the cluster interaction 
by interaction using the rejection method: we select two 
black holes randomly, compute the cross section E of the 
interaction, divide by the largest possible cross section 
E max (which is the cross section of capture by the two 
most massive black holes in the sample), and then com- 
pare that ratio with a uniform random deviate x € [0, 1). 
If x < E/E max we accept the interaction, otherwise we 
draw again. 

If the intera ction is between two black holes, then we 
use the recent lLousto et alj (|2012l ) formula for the kick. 
If the kick is greater than the escape speed v csc = 4<j 
(typical of a core-collapsed cluster) we assume that the 
remnant has been ejected from the cluster and thus we 
remove both black holes from the sample. Otherwise, we 
assume the remnant remains, hence we sum the masses of 
the holes and estimate t he spin of the remnant following 
the prescription given in iRezzolla et ail ([20081 ) . 

Figure [2] shows the results. Here we plot the fraction 
of clusters that retain a black hole after subcluster evo- 
lution, and the median mass of the final black hole if 
one remains, as a function of the velocity dispersion a of 
the cluster. For each velocity dispersion we performed 
10 4 simulations. For low escape speeds, almost all merg- 
ers between black holes eject the remnant, hence reten- 
tion depends on whether the initial number of holes is 
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even or odd. As the escape speed increases, so does the 
probability that a merger will not eject the remnant; for 
w C sc SS 100 km s _1 it is most probable that this happens 
when the spins of the holes are low and their masses are 
close to each other (note from symmetry that there is 
zero recoil from the merger of equal-mass nonspinning 
holes). As the escape speed increases further, mergers 
between black holes of different masses can be retained, 
until at v csc ^> 800 km s -1 a runaway occurs and a single 
victorious black hole is usually the result. 




Velocity Dispersion (km s" 1 ) 



Fig. 2. — Fraction of clusters of a given velocity dispersion that 
retain a black hole after a succession of mergers (upper left curves) 
and, if a black hole is left, the median mass of the remaining black 
hole (lower right curves). For this figure we ignore the effects of 
hard binaries formed by the interactions of three initially hyper- 
bolic black holes (see text), hence there is a difference between 
cases with an initially even and an initially odd number of black 
holes. The solid curves are for 100 initial black holes, and the dot- 
ted curves are for 101 initial black holes; the asymmetry in retained 
fraction at low velocity dispersions is because if every black hole 
merger results in an ejection, an initially even number will leave 
behind no black holes whereas an initially odd number will leave 
behind one. We assume an escape speed that is four times the 
velocity dispersion. This figure demonstrates that retained run- 
away mergers leading to massive seeds are only likely for velocity 
dispersions Jj 200 km s — . 



From these simulations we can argue that for clusters 
with velocity dispersions <J 100 km s _1 a runaway is un- 
likely, but that there is roughly an equal chance of leav- 
ing behind either one or zero holes (depending largely on 
the parity of the initial number until a > 60 km s _1 ). 
When there is a black hole left behind it is likely for 
a < 100 km s _1 to be at the low end of the mass distri- 
bution (~ 5 Mq), because such black holes are initially 
more common. In addition, lower-mass black holes have 
a lower cross section for capture and hence an enhanced 
probability of survival. 
The subsequent evolution has two possibilities: 
One black hole remains. — Then as we discuss in the 
next section, the black hole will sit near the center of the 
high number density distribution of stars. Tidal disrup- 
tions will add a few tens of percent of the stellar mass to 
the hole, mostly within a few weeks or less of the initial 
disruption, hence the hole will grow quickly. Given that 



interactions with the stars cannot eject the hole from the 
cluster, it will become a massive black hole in a short 
timescale. 

No black holes remain. — In this case, the stars will un- 
dergo runaway collisions with themselves, leading to the 
production of a massive star that will then beco me a 
black hole (e.g.. lPortegies Zwart fc McMillan|[2002t ) . The 
situation then reduces to the previous case, because the 
time needed to produce a second black hole, which could 
potentially eject the first, is significantly larger than the 
time needed for the first hole to increase its mass to the 
point that it can no longer be ejected. 

We now discuss these possibilities in greater depth. 

3.1. Interactions between stars and a black hole 

Although the central density after core collapse is for- 
mally infinite, the finite number of stars means that this 
translates to a few stars in a small region near the core. 
For example, if we consider the inner ~ 10 solar-type 
stars after core collapse and continue to assume a con- 
stant velocity dispersion, then they are in a region r = 
GM/a 2 ~ 5 AU(M/10 M )(ct/4O km s" 1 )" 2 in radius, 
with a resulting number density of n > 10 14 pc -3 . Even 
the inner 1000 stars are in a region with n > 10 10 pc~ 3 , 
so interactions will be common and rapid. 

Stellar tidal disruption by black holes. — A promising 
mechanism for such runaway growth is tidal disruption of 
stars by stellar-mass black holes. The critical pericenter 
for tidal disruption of a star of mass m and radius R by 
a black hole of mass M is 

rp.tidal = (SM/my^R . (11) 

Thus the gravitationally focused cross section for tidal 
disruption, assuming that the black hole mass greatly 
exceeds the stellar mass, is 

Stidai ~ 10 26 cm 2 (M/10 A/ Q ) 4/3 (cr/40 km s" 1 )" 2 (12) 

for solar-type stars. This is roughly an order of mag- 
nitude greater than the star-star collision cross section 
discussed later, and two orders of magnitude larger than 
the black hole - black hole capture cross section. More- 
over, the rate is nonlinear in the mass of the black hole 
(E oc Af 4 / 3 ). Thus the conditions for a runaway exist. 

If tidal disruption does occur, then the mass will 
be force fed to the black hole at an extremely super- 
Eddington rate. Studies suggest that fallback initially 
occurs over several ti mes the internal dynamic al time 
of the disrupted star (|Evans fc K ochanckl ll989l ). which 
is several hours for a solar-type star. The accretion 
rate is therefore many millions of times the Edding- 
ton rate of a stellar-mass black hole. Analyses of 
such supercritical accretion (e.g.. IMaraschi et al.l 119761: 
Begelmanl 119791: [Jaroszvnski et alj Il980t IPopham et al.l 
1999; lOhsuea et al. 2005) indicate that the matter will 
indeed flow into the hole at that rate, but that most of 
the photon luminosity that is generated will be advected 
in with the very optically thick matter (hence although 
the accretion rate is tremendously super-Eddington, the 
luminosity could be limited to Eddington or slightly 
higher). Thus it is expected that within a matter of days, 
i.e., much shorter than any other relevant timescale, most 
of the bound remainder of the star will flow onto the hole. 
If this is the case, then the majority of the accretion will 
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finish without harassment from additional encounters by 
stars. If, on the contrary, the accretion rate is actually 
limited to the Eddington rate then the time needed to 
accrete most of the matter is much longer than the time 
to the next encounter, and the disk might be disrupted, 
leading to negligible growth of the hole. 

The unbound remnant of the star will be thrown out- 
wards at speeds comparable to the orbital speed at tidal 
disruption, which is ~ 800 km/s(M/10 Mo) 1 / 3 for a 
solar-type star. This is much larger than the escape 
speed, so the wind will depart ballistically unless it runs 
into many times its own mass in gas in the cluster. How- 
ever, given that the virial temperature of the cluster is 
~ 10 5 K(cr/40 km s" 1 ) 2 and that cooling is extremely 
efficient at that temperature, the total amount of gas in 
the cluster at a given time will be small even though its 
escape speed is sufficient to retain winds from red giants 
or (earlier, when more massive stars existed) planetary 
nebulae. Thus we assume that the unbound gas simply 
escapes from the cluster. The ratio of unbound gas to 
gas that accretes onto the black hole is rather uncertain. 
The initial disruption leay es about half the mass bound 
(E vans fc Kochanekl 119891) . but shocks upon the return 
of the bound mat ter might unbind additiona l mass. In 
a recent study by iStrubbe fc Quataertl (|2009f ) they con- 
sider different ejection fractions ranging from / osc = 0.5 
(corresponding to negligible return shocks) to / csc = 0.8 
(corresponding to powerful return shocks). In our sce- 
nario, the upshot is that because a single black hole will 
grow, its growth will eject up to a few times its own mass 
in stellar debris. Until this reaches at least hundreds, and 
probably thousands, of solar masses this will be such a 
small fraction of even the core mass that we expect it to 
have a minor effect on the dynamical evolution. 

Star-star collisions. — At the velocity dispersions we 
consider, these collisions are likely to lead to merg- 
ers with little mass loss, because a ~ 40 km s _1 is 
much less than the escape speed ~ 600 km s _1 of 
a solar-type star. For the same reason, these col- 
lisions are gravitationally focused, with a cross sec- 
tion E = 7rr 2 (l + 2GM tot /(r p a 2 )) » 2ir(GM tot /a 2 )r p 
for a pericenter distance r p and a total mass be- 
tween the stars of M to t- The relevant pericenter dis- 
tance is the sum of the stellar radii, which is 2Rq ps 
0.01 AU for two solar-type stars, hence for two such 
stars E » 1.5 X 10 25 cm 2 (cr/40 km s^ 1 )" 2 . The char- 
acteristic time of interaction is then r — l/(nE<r) « 
10 6 yr(n/10 10 pc~ 3 )- 1 (M tot /2 Mq)- 1 ^/^ km a -1 )" 1 . 
Note that as a result even for the inner ~ 10 3 stars the 
collision time for solar-type stars is much less than their 
~3x 10 7 yr Kelvin-Helmholtz time, hence the stars will 
not be able to radiate their collisional energy before the 
next collision. However, because the velocity dispersion 
is < 0.1 times the stellar escape speed, the energy added 
is minor and most of the pressure holding up the colli- 
sion product stems from gravitational contraction rather 
than either collision energy or nuclear energy; these are 
thus not stars in the standard sense, and need not have 
luminosities as high as those of main sequence stars of 
the same mass. 

In addition, on the main sequence, stellar radii increase 
with increasing mass, hence the rate of interactions in- 



creases more than linearly with increasing stellar mass. 
An additional factor is that more massive stars tend to 
sit closer to the center of the potential, where the num- 
ber density of objects is greater. The conditions are thus 
ripe for a runaway, and indeed runaway merging of stars 
has been proposed as a mechanism for the generation of 
supermassive stars that later evolve into interm e diate- 
mass black holes ()Portegies Zwart fc McMillan! I2002t 
iFreitag et al.ll2006l) . It has been suggested that the high 
wind rates expected for high-mass stars can severely limi t 
the growth of supermassive stars (jGlebbeek et al.l l2009) . 
Note, however, that these wind rates are based on ex- 
trapolations of winds for main sequence stars, and as in- 
dicated above the collision products will be substantially 
larger and less luminous than main sequence stars. In- 
deed, the collision products are more likely to be a "bag 
of cores" than an actual star, where an extended gaseous 

envelope engulfs an ensemb le of stellar core s. 

We do note that although iGlebbeek et all (|2009f) argue 
that winds may prevent the formation of intermediate- 
mass black holes, they find that runaway collisions pro- 
duce stars massive enough to evolve to normal stellar- 
mass black holes, at least. Thus for our purposes we 
assume that star-star collisions will lead to black hole 
production. 

The question is then whether the first black hole that 
forms has enough time to consume many stars so that 
by the time the next black hole forms, the first one is 
so massive that any BH-BH merger will produce a weak 
recoil that retains the remnant in the cluster. We argue 
that this is in fact the case: the first black hole to form 
will be at the center of the mass distribution, where the 
number density is the highest. If this is, for example, in 
the region occupied by the inner ~ 100 stars, then the 
number density is such that tidal disruptions of stars by 
even a 10 M@ black hole occur on average once per few 
hundred years, and the interaction time scales as M~ 4 / 3 . 
Thus the hole will double its mass every few thousand 
years, i.e., in a time vastly shorter than the lifetime of 
even the most massive stars. We note that although the 
segregation of the black holes to the center of the cluster 
and their ejection leads to some flattening of the stellar 
number densi ty near the center of the cluster (see, e.g., 
IAarsethll2012l and in particular his Figure 8 for a recent 
N-body simulation), we expect that when most of the 
black holes have been ejected the stars near the core, 
which have a very short relaxation time, will regrow the 
cusp. Hence the first black hole formed due to runaway 
stellar collisions will be able to increase its mass by a 
large factor before any other new-generation black hole 
forms. 

3.2. Minimum mass of central BH 

We consider here the evolution of a cluster with a ve- 
locity dispersion large enough to guarantee core collapse. 
If a black hole grows in the cluster, what is a rough ap- 
proximation to its minimum mass? We will approach 
this question in two different ways. First, we will de- 
termine the mass of a black hole nailed to the center of 
an n oc r~ 2 core collapse cluster such that dynamical 
processes around the black hole can supply enough heat 
to help forestall further core collapse. Second, we will 
apply the criterion that the wander radius of the black 
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hole must be less than its radius of influence, under the 
assumption that otherwise the number of stars bound 
to the black hole would be much less, hence its heating 
influence would be reduced. 

We will assume as before that the mass of the nu- 
clear star cluster is related to its velocity dispersion by 
M c i ~ 10 6 M (ct/4O kms -1 ) 4 . For a core collapse 



mass (from energy equipartition arguments), hence for 
this system we expect 



cluster with 



', the velocity dispersion is the 



same at all radii and the gravitational binding energy 
is -Ebind = (l/2)Aicicr 2 from the virial theorem. 

This energy must be compared with the available en- 
ergy (as defined before) from dynamics around the cen- 
tral black hole. The available energy per unit mass 
around the black hole that we found previously is 
GM-BB,/(^ r p,min)t where r Pimm is the minimum pericen- 
ter distance of an orbit that can last long enough for 
significant dynamical interactions. For a black hole, 
the relevant time is the time for gravitational radia- 
tion to cause the object to spiral in; this time scales 
as T ~ (mM| H ) r 4 , roughly, so for a fixed T we have 

r p ,mm oc Af^ 2 . We used r p , min w 0.01 AU for 10 M Q 
(giving an inspiral time of a few million years), so we 
will adopt r p , min = 0.1 AU(M BH /10 3 Mo) 1 / 2 . 

If the distribution of stars around the black hole is a 
steep cusp, then the stellar mass in the radius of influence 
of the black hole equals the mass of the black hole (this 
need not be true if the den sity distribution ha s a core 
profile; see equation (14) of lLasota et alJl201.ll) . When 
we compare the available dynamical energy of the stars 
around the black hole with the binding energy of the 
cluster, we find that 



Mbh £ 500 M Q (cr/40 km s" 1 ) 



(13) 

is required for the stars around the black hole to provide 
sufficient energy to hold off collapse. 

We can also approach this from a different angle. A 
finite-mass black hole will not be nailed to the center 
of the cluster. Instead, it will wander due to stochastic 
dynamical interactions. If the wander radius is less than 
its radius of influence then we can suppose that it is near 
the center of the stellar distribution where encounters 
are frequent, but if the wander radius is larger then this 
need not be the case and heating could be less efficient. 
Thus a different criterion is r wan d CT < Hnfl.BH- Suppose 
that there is a nearly constant-density core in the inner 
10% of the cluster; then the scale height of a species in a 
cluster is inversely proportional to the square root of its 



^wander 



r cl ((m)/M BH ) 1/2 



(14) 



Here (m) is the average mass of a star. The cluster radius 
is r c \ = GM c \/a 2 and the radius of influence of the black 
hole is GMbh/c 2 , hence the wander criterion is 

a/2 



0.1(GM cl /a 2 )((m)/M BH y 
Mbh 
M bh 



£ GMbh/ct 2 
£ (0.01 M^(m)) 1 / 3 
;> 2 x 10 3 M q (ct/40 kms" 1 ' 
(15) 

where in the last line we assume (m) = 1 Mq. 

Recall that these are lower limits on the mass of the 
central black hole. The mass could be considerably 
greater depending on long-term accretion of stars or gas. 

4. DISCUSSION AND CONCLUSIONS 

We have discussed the evolution of a relaxed cluster 
that has a velocity dispersion a ^> 40 km s -1 , which is 
large enough to render binaries insignificant, but that 
does not initially contain a massive central black hole. 
We argue that a massive hole will inevitably form if it 
can swallow tidal debris rapidly: interactions in the black 
hole subcluster will leave either zero or one hole. In the 
case of zero, a black hole will form from the product of 
runaway stellar merging. In either case, the hole will 
feed quickly from the remnants of the stars it tidally dis- 
rupts, and hence will grow until it has significant dynam- 
ical effects on the cluster and thus slows its own growth. 
It is not guaranteed that the holes will then follow the 
same M — a relation that exists for higher velocity dis- 
persion systems. It is also not guaranteed that clusters 
with lower velocity dispersions will not have black holes, 
but it is possible that massive black-hole formation is 
prevented as long as bina ries have a significant heating 
effect (see lGill et al J 120081 for a numerical exploration of 
the heating due to binaries or a massive central object). 
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